Universal charge-radius relation for subatomic and astrophysical compact objects 
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Electron-positron pair creation in supercritical electric fields limits the net charge of any static, 
spherical object, such as superheavy nuclei, strangelets, and Q-balls, or compact stars like neutron 
stars, quark stars, and black holes. For radii between 4 x 10 2 fm and 10 4 fm the upper bound on 
the net charge is given by the universal relation Z — 0.71i?f m , and for larger radii (measured in fm 
or km) Z = 7 x 10~ 5 i? 2 m = 7 x 10 31 i?^ m . For objects with nuclear density the relation corresponds 
to Z w 0.7A 1/3 (10 s < A < 10 12 ) and Z « 7 x 10~ 5 yl 2/3 (A > 10 12 ), where A is the baryon number. 
For some systems this universal upper bound improves existing charge limits in the literature. 

PACS numbers: 04.70.-s, 21.65.Mn, 21.65. Qr, 97.60.Jd 



Spontaneous formation of real electron-positron pairs 
in supercritical electric fields is known to lead to screen- 
ing of highly charged objects [j], Q. This is in addition 
to the vacuum polarization effect caused by virtual pairs 
In the following it will be shown, that pair forma- 
tion implies a universal upper limit on the net charge 
of any static, spherical object of given radius. The up- 
per bound on the net charge for R > 4 x 10 2 fm is 
Zoo = 0.71-Rfm = 7.1 x 10 17 i?km, where the radius is 
measured in fm or km to show the typical scales for 
subatomic objects or compact stars like neutron stars, 
quark stars and black holes, and the subscript oo indi- 
cates infinite time available for pair formation. Hypo- 
thetical superheavy nuclei, quark nuggets (also known as 
strangelets), and other objects with nuclear matter den- 
sity have -Rf m « A 1 / 3 , where A is the baryon number, so 
the charge bound can be recast as Zoo ~ 0.7 A 1 / 3 . This 
universal charge bound is complementary to and often 
significantly more restrictive than other bounds in the 
literature. However, it implicitly assumes that infinite 
time is available to populate electron levels via pair cre- 
ation. Taking the relevant timescales into account the 
maximum charge increases for R > 10 4 fm (A > 10 12 ) to 
Z « 7 x KT 5 i? 2 m « 7 x lCr 5 ^ 2 / 3 (the latter expression 
assuming nuclear matter density). This relation improves 
existing limits for e.g. strangelets, reproduces an earlier 
result for massive black holes, but is inferior to charge 
limits based on stellar stability for gravitationally bound 
neutron stars when only finite time is available for pair 
creation. 

In the following the upper bound on the charge of 
spherical objects will be derived in the context of a rel- 
ativistic Thomas-Fermi model calculation following the 
approach of Miiller and Rafelski [l]. The universal re- 
lation will be analytically derived in the high mass, 
infinite time limit, and as typical examples homoge- 
neously charged superheavy nuclei and color-flavor locked 
strangelets dominated by surface charge will be shown 
numerically to approach the universal charge-radius re- 
lation. Finite time for pair creation will then be taken 



into account and the corresponding (higher) charge limits 
derived. The results will be compared to earlier charge- 
mass relations for strangelets, and to limits for astro- 
physical objects. Whereas the derivations do not include 
general relativity it is argued, that the universal charge- 
radius limit should also apply approximatively to black 
holes, and in fact the results agree with a limit derived 
for black holes in general relativity. 

In a continuous approximation at zero temperature 
(more about finite temperature effects later) the number 
density of electrons, n c , is given by the electron mass, 
rn , and Fermi energy, /z , via 
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Here the effective chemical potential /x° ff = fj. e — eq 
where <j> is the electric potential, and the ^-function takes 
into account that /x e must exceed the electron mass. In 
the usual Thomas-Fermi model describing neutral atomic 
systems one takes the energy needed to add an additional 
electron, /xjj ff = m e . In the present context the focus is on 
maximally charged rather than neutral systems, and here 
one instead takes /ig ff = — m e , corresponding to the top 
of the negative energy sea [jij . This represents a situation 
where all levels accessible to spontaneous vacuum decay 
are filled, and as shown in 0, Q it reproduces the results 
of single-particle calculations carried out for core charges 
up to a few hundred that demonstrate how more and 
more real electron states dive into the negative energy 
continuum as the core charge increases, leading to the 
creation of a negatively charged vacuum. Thus 
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In regions of space containing only electrons Poisson's 
equation is given by V 2 (e<^)) = 47re 2 n c , or in dimension- 
less units, where the radial coordinate £ = m e r, and the 
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normalized potential y = e<p/m c , with a 



V 2 y = ^ [y 2 - 2y] 
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In regions with additional charge, a corresponding source 
term must be added to the right-hand side of Poisson's 
equation. Two illustrative types of core charge distribu- 
tions have been considered: A uniformly charged spheri- 
cal core (like a nucleus) and a uniformly charged spheri- 
cal shell with no internal charge (like a color-flavor locked 
quark matter lump or a large, perfectly conducting nu- 
clear matter system, where the net core charge moves to 
the surface). Together these idealized distributions span 
the likely range of real distributions, and as shown in the 
following, they lead to identical results for the screened 
charge in the limit of large system radius. 

Boundary conditions for Poisson's equation are Vy — > 
for £ — > and y — > for £ — > oo. But the actual behav- 
ior of y for large £ is explicitly known in the maximally 
charged case. Denote the total (positive) core charge by 
Z COIC and the total number of electrons by N e . Then 
Zoo = Zcote — N e > is the net charge of the maximally 
charged system seen by an observer outside the radius 
£2, where the electron density drops to zero, 1/(^2) = 2 
(corresponding to e0 = 2m e ). But to this observer 
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Using this relation for £ = £ 2 gives 

2m c r2 



(4) 
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where r 2 denotes the physical radius corresponding to 
£ = £2- In the bulk of a homogeneously charged core 
electrons neutralize the local core charge, and a devia- 
tion from local charge neutrality occurs only very close to 
the surface. As confirmed by numerical solutions of Pois- 
son's equation the characteristic width of the screening 
electron cloud outside the core charge is of order m^ 1 or 
a few hundred fm (deviations from local charge neutral- 
ity occurs within a similar zone inside the surface) . Thus 
in the limit of r 2 S> m" 1 the width is small compared to 
the physical radius of the core, R, and to good approxi- 
mation the net charge of a maximally charged spherical 
system with R 3> m^ 1 w 4 x 10 2 fm is therefore (since 
R » r 2 ) 
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= 0.71R{ m = 7.1 x 10 17 i? k 
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For objects with density of order nuclear matter den- 
sity i?f m w A 1 / 3 . Therefore, for such objects charge and 
baryon number are related by 



Figure 1 shows how well the universal charge-radius re- 
lation fits numerical solutions to the relativistic Thomas- 
Fermi model at large radius. The example here is for 
color-flavor locked strangelets composed of a core of equal 
numbers of up, down, and strange quarks, leading to 
zero quark charge density in the bulk of the core, but 
with a core surface charge due to surface depletion of 
strange quarks, so that Z corc ps 0.3A 2 / 3 [3j and radius 
Rim ~ 1.1 A 1 ' 3 (such a system has a discontinuity in 
dy/d£ at the core surface in addition to the boundary 
conditions previously described for the Poisson equa- 
tion). One clearly sees how the actual solution of Pois- 
son's equation for the net screened charge shifts from 
Zoo ~ Z COTC to the universal relation Z m « 2£ /a when 
going from small (£ -C 1) to large (£ » 1) radii. Sim- 
ilar behavior is found for systems with other Z COIC (A) 
dependence and regardless of whether the core charge is 
uniformly distributed, or distributed as a surface charge 
as in the case of strangelets and for cores that behave 
as ideal conductors. In all cases tested, the maximal 
charge approaches the same universal charge-radius re- 
lation Zoo ~ 0.71iifm for m c R = £ 3> 1 (physical radius 
R exceeding 4 x 10 2 fm) as expected. 

So far the calculations have assumed a static situation 
with infinite time available to fill all accessible electron 
states. However the pair formation process involves tun- 
neling and the rate for this (number of pairs produced 
per volume per time) has been shown to be Q 
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where the critical field is E cr — m^/a 1 / 2 ~ 1.3 x 
10 16 V/cm. Taking E = Za 1 / 2 / R 2 to represent the "sur- 
face region" of the spherical charge distribution studied 
here, and including only the leading n = 1 term in the 
sum, this corresponds to 
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Zoo ~ 0.7 A 1/3 for A > 10 s 



(7) 



Taking the active charge producing layer to have a vol- 
ume V = 4nR 2 AR 1 with thickness AR * m e 1 , assuming 
WV = —dZ/dt, and defining a characteristic timescale r 
for charge equilibration from dZ/Z = —dt/r, this leads 
to a timescale 

t = —^—xexp(x) — 5.5 x 10 _19 sxexp(a;), (10) 

where x = 7r£ 2 /aZ. Because of the exponential a fi- 
nite time available for pair production corresponds to a 
(roughly) fixed value of x for a wide range of t (it also 
means that the results do not depend crucially on the 
assumptions about V; taking V = AnR 3 /3 instead of the 
surface layer chosen makes little difference to the charge 
relations below). In other words, pair production in a fi- 
nite time gives a maximum charge value which is propor- 
tional to £ 2 , rather than £ as were the case with infinite 
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time available. For example, for r = 1 s the maximally 
allowed charge becomes 

Z ls « 11. 2£ 2 « 7 x W~ 5 R{ m ~7x 10~ 5 A 2/3 , (11) 

where the last equality assumes nuclear matter den- 
sity. For r = 10 10 years one gets very similar re- 
sults, except for dividing the numerical prefactors by 2.0, 
whereas a typical weak interaction timescalc of 10~ 10 
s corresponds to multiplication of the numerical pref- 
actors by 2.4. Concentrating on Z\ s we therefore have 
two regimes for the universal charge-radius relation. For 
Z\ S {C) < -^oo(^) pair production is sufficiently rapid 
to fill the electron states as assumed in the relativistic 
Thomas- Fermi model, and therefore Z « Z^ w 2£/a for 
1 < £ < 25. For Z la (g) > ^oo(C) (6 > 25) we have instead 
Z m Zi s ~ 11.2£ 2 . The corresponding division between 
the two regimes in terms of radius and baryon number 
(assuming nuclear matter density) is Rf m ~ 1.0 x 10 4 
and A w 10 12 respectively. The higher charge permit- 
ted for macroscopic spheres because of the finite time 
effect explains why a Van de Graaff generator can work 
at potential in excess of 1 MV as would be the limit with 
infinite time to screen the charge with electrons from pair 
production. 

The existence of a universal maximal charge-radius re- 
lation for spherical objects regardless of their physical 
nature is interesting in its own right, but it also has ap- 
plications to several areas of subatomic physics and as- 
trophysics as outlined in the following. 

It has been speculated that there could exist branches 
of mctastable superheavy nuclei. If so, the charge of the 
maximally ionized superheavy ion should obey the rela- 
tion Z w 0.7A 1 / 3 (10 8 < A < 10 12 ) or Z « 7 x 10- 5 A 2 / 3 
(A > 10 12 ). Because the nuclear charge increases almost 
linearly with A, charge screening becomes important al- 
ready for A « 10 3 (as also pointed out in [l[), but numer- 
ical solutions of Poisson's equation show that one needs 
£ > 1 or A > 10 8 for the universal relation to become 
quantitatively accurate. 

The possible existence of metastable or even stable 
quark nuggets or strangelets has been widely discussed 
[5j. For large chunks of quark matter, which could exist 
in our Galaxy as a result of binary compact star colli- 
sions, the (Z, A)-relation in the high-mass limit has so 
far been taken as either Z w 8A 1 / 3 (for non-color-flavor 
locked strangelets with A < 10 7 [6]), or Z « 0.3A 2 / 3 (for 
color-flavor locked strangelets regardless of A Q). The 
possible importance of charge screening due to supercrit- 
ical field electrons was mentioned by Farhi and Jaffe Q, 
and a first numerical study was included in 8]. From 
the discussion above it follows that the universal relation 
Z « 0.7A 1 / 3 (10 8 < A < 10 12 ) and Z « 7 x 10~ 5 A 2 / 3 
(A > 10 12 ) can be applied for A > 10 8 as an upper enve- 
lope on the strangelet charge regardless of the details of 
the quark phase, and for color-flavor locked quark matter 
this envelope represents the actual maximum ionization 



state possible, since the "envelope" value of Z is smaller 
than the core charge 0.3A 2 / 3 . 

Q-balls have been suggested in various varieties, some 
of which can be charged [9j . Again, the universal relation 
derived here should apply. 

Turning to astrophysical objects, it is normally as- 
sumed that compact stars such as neutron stars and 
quark stars are close to electrical neutrality, since a net 
positive charge would attract electrons from the inter- 
stellar medium. However, many papers have dealt with 
the theoretical possibility of non-zero charge and its in- 
fluence on stellar structure, and limits (typically from ar- 
guments related to stellar stability) have been placed on 
the charge allowed. From a comparison of gravitational 
binding and electric repulsion, several authors (e.g. (lol |) 
have found that the net charge of gravitationally bound 
stars is limited by Z < 10~ 36 A, where the baryon num- 
ber A w 10 57 for a typical compact star, thus Z < 10 21 . 
From the universal relation derived above, we have the 
limit Zis < 7 x 10~ 5 A 2 / 3 , which is conceptually very dif- 
ferent, but happens to give a larger number for A ss 10 57 , 
namely Z± s < 7 x 10 33 . This shows that stabilization 
of charged neutron stars should result from trapping of 
real electrons from the surroundings; there is not enough 
time to ensure sufficient numbers of electrons from pair 
creation. With infinite time available, the static maxi- 
mal charge would be Z^ w 0.7A 1/3 w 10 19 , i.e. stable 
with respect to the electric repulsion, but to approach 
this limit requires orders of magnitude more time than 
the age of the Universe if one were to start out with a 
much higher charge. 

Strange stars consisting of quark matter are self-bound 
due to strong interactions in addition to gravity, so the 
relevant stability limit here comes from a comparison of 
the Coulomb energy and the total (strong interaction) 
binding energy per baryon. Stability requires Z 2 a/AR 
to be less than a few MeV, which is easily satisfied (by 8 
orders of magnitude) for Z = Zi s . Therefore the charge 
of strange stars may in principle saturate the universal 
limit Zis ~ 7 x 10 -5 A 2 / 3 , and such a system could have 
interesting properties related to pair creation even at zero 
temperature in addition to the (much faster) finite tem- 
perature pair creation phenomena discussed in [llj ]. 

Black holes are formally exempt from the treatment 
above since they are described by general relativity. How- 
ever, formation of an astrophysical black hole involves the 
collapse of a mass (and for the present purpose charge) 
that seen from an outside observer only asymptotically 
reaches the horizon. Therefore one might expect that the 
universal charge-radius relation would indeed be obeyed 
by the black hole because it would be obeyed during the 
formation process, and the maximal net charge seen from 
the outside again should not exceed 7 x 10 31 i? 2 m . Stud- 
ies of pair formation in a general relativistic description 
of black holes has in fact led to results almost identical 
to this [l2T |. a result which is many orders of magnitude 
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below other limits derived for the maximal charge of a 
stable black hole as long as the mass is below 10 8 solar 
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masses 

The derivations above were all based on zero temper- 
ature relations. These relations remain valid for temper- 
ature T <C fJ-e ~ m c = e<j) — 2m c , so at low tempera- 
tures thermal electrons can be neglected except near the 
edge of the zero temperature maximally charged configu- 
ration, where e<j) — > 2m c . The lowest order thermal con- 
tribution to the electron number density in this regime is 
roug hly ml /2 T^ 2 , giving a thermal contribution to the 
total charge which is suppressed by (T/m ) 3 / 2 relative 
to the zero temperature charge Z± s . Therefore the zero 
temperature relations provide good approximations for 
T <c m e . 

Another caveat in applying the universal charge-radius 
relation derived above to real subatomic or astrophysi- 
cal objects is related to the tunneling timescale involved. 
The larger the system, the lower is the probability of fill- 
ing all electron levels on a short timescale. Therefore, 
even though the universal relation is independent of the 
internal composition and structure of the spherical core, 
the realization of the universal behavior in realistic sys- 
tems will depend somewhat on the nature of the object 
and its previous history. Systems with a high electri- 
cal conductivity (e.g. nuclear matter or quark matter) 
will rearrange any initial net charge in a way so that the 
charge gets concentrated in a thin surface layer. The 
small width of this layer and the surrounding electron 
atmosphere improves chances for the universal charge- 
radius relation to be realized. Other systems, such as 
strangelets or quark stars composed of color-flavor locked 
quark matter, are electrically neutral in the bulk quark 
phase already 11 311 - and have only a net quark charge on 
the surface [3, LL4j , so here the electron shield is thin al- 
most "by construction" . Positively charged objects in 
general tend to neutralize by trapping electrons from the 
surroundings, e.g. from the interstellar medium. There- 
fore, most of the bulk of macroscopic objects will be close 
to charge neutral, and the interesting physics will be re- 
lated to processes ionizing such objects. Such ionization 
will take place "from the outside" and involve a restricted 
radial range, again minimizing the timescale problem. 

In conclusion, it has been demonstrated that electrons 
created in supercritical fields lead to a universal relation 
between the maximal charge and radius of any static, 
spherically symmetric object with a size exceeding a few 
hundred fm. The relation limits the charges of objects 
such as superheavy nuclei, strangelets, Q-balls, neutron 
stars, quark stars, and black holes. For some of the ob- 
jects studied, the new limits are more restrictive than 
other charge relations and limits existing in the litera- 
ture. 
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FIG. 1: Charge as a function of scaled radius, f. The lower 
(blue, dotted) line represents the universal charge-radius re- 
lation Zoo ~ 2£/a valid for large radii and infinite time. The 
upper (green, dashed) line is the quark charge of color-flavor 
locked strangelets, and the middle points with fitted (red) 
curve shows the maximal net charge, Z, as calculated nu- 
merically from the relativistic Thomas-Fermi model. As can 
be clearly seen, the electron screening due to the supercritical 
electric field becomes important for £ > 0.1, and the universal 
charge-radius relation is a good representation for f 3> 1. 



